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Context. The distance of the very young open cluster Westerlund 2, which contains the very massive binary system WR 20a and is 
likely associated with a TeV source, has been the subject of much debate. 

Aims. We attempt a joint analysis of spectroscopic and photometric data of eclipsing binaries in the cluster to constrain its distance. 
Methods. A sample of 15 stars, including three eclipsing binaries (MSP 44, MSP 96, and MSP 223) was monitored with the FLAMES 
multi-object spectrograph. The spectroscopic data are analysed together with existing B V photometry. 

Results. The analysis of the three eclipsing binaries clearly supports the larger values of the distance, around 8 kpc, and rules out 
values of about 2.4 - 2.8 kpc that have been suggested in the literature. Furthermore, our spectroscopic monitoring reveals no clear 
signature of binarity with periods shorter than 50 days in either the WN6ha star WR20b, the early O-type stars MSP 18, MSP 171, 
MSP 182, MSP 183, MSP 199, and MSP 203, or three previously unknown mid O-type stars. The only newly identified candidate 
binary system is MSP 167. The absence of a binary signature is especially surprising for WR 20b and MSP 18, which were previously 
found to be bright X-ray sources. 

Conclusions. The distance of Westerlund 2 is confirmed to be around 8 kpc as previously suggested based on the spectrophotometry 
of its population of O-type stars and the analysis of the light curve of WR 20a. Our results suggest that short-period binary systems 
are not likely to be common, at least not among the population of O-type stars in the cluster. 
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1. Introduction 

The open cluster Westerlund 2 lies in a blowout region of the 
giant Hn complex RCW49. Interest in the stellar population 
of Westerlund 2 was triggered by two independent observa- 
tional studies. On the one hand, RCW49 was observed with the 
Infrared Array Camera aboard Spitzer in the framework of the 
GLIMPSE legacy survey (Churchwell et al. 120041 . These obser- 
vations revealed strong evidence of ongoing star formation activ- 
ity and underlined the need for a detailed study of the properties 
of the early-type stars in the cluster core and their impact on the 
surrounding nebula. On the other hand, WR 20a - one of the two 
Wolf-Rayet stars in RCW49 - was found to be a very massive 
eclipsing binary consisting of two WN6ha stars with individual 
masse s of about 80 M G (Rauw et al. 120041 12005 ; Bonanos et al. 
2004). Since previous spectroscopic studies revealed only one 
06: and six 07: stars, in the cluster core (Moffat, Shara & Potter 
119911 ). the existence of a pair of stars that massive was somewhat 
unexpected. 

The nature of Westerlund 2 thus remains to be established: 
is it a super-cluster maybe harbouring dozens of very massive 
early O-type stars, or is it "just" a medium-size agglomerate 
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that produced the very massive WN6ha + WN6ha binary sys- 
tem WR 20a more or less by chance?[^This question called for a 
re-investigation of the population of early-type stars in this clus- 
ter. A first step in this direction was performed by Rauw et al. 
(120071 hereafter Paper I), who presented a photometric monitor- 
ing of the cluster as well as a spectroscopic snapshot study of 
its brightest members. This study allowed us to discover three 
previously unknown eclipsing binaries among the likely clus- 
ter members and to derive spectral types between 03 and 06.5 
(i.e. significantly earlier than previously thought) for the twelve 
brightest O-stars in the cluster. Moreover, we inferred a spec- 
trophotometric distance of 8.0± 1 .4 kpc. While this distance mea- 
surement is in excellent agreement with that of WR 20a (Rauw et 
al. 120051 ). it is nevertheless problematic in several ways. Indeed, 
despite the much earlier spectral types inferred for the most mas- 
sive cluster members, the current census of early-type stars in 
Westerlund 2 accounts for only 20% of the ionizing photons re- 
quired to explain the observed radio emission of the RCW49 
complex. Moreover, a near-IR photometric study of the clus- 
ter by Ascenso et al. (120071 ) resulted in a distance estimate of 
2.8 kpc, inferred from a comparison between their JHK S data 
with pre-main sequence evolutionary tracks. The origin of these 
discrepancies remains unclear. More recently, yet another dis- 
tance determination was obtained. From the mean velocity and 



1 Indeed, the existence of a relation between the mass of an open clus- 
ter and that of its most massive member remains disputed (see Weidner, 
Kroupa & Bonnell 2010 and references therein). 
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velocity spread of molecular gas in the molecular cloud associ- 
ated with RCW 49, Furukawa et al. d2009) inferred a kinematic 
distance of 5 . 4+ J ^ kpc. These authors further suggested that the 
collision between two molecular clouds might have triggered the 
formation of the Westerlund 2 cluster. 

Determining the correct distance of this cluster is a fun- 
damental issue: if the distance were confirmed to be 8.0 kpc, 
Westerlund 2 might be one of the most massive clusters of our 
Galaxy (along with NGC3603, CygOB2, and Westerlund 1). 
Otherwise, if the distance were established to be much shorter 
than 8.0 kpc, WR20a would be left in relative isolation, thereby 
challenging the theories of competitive accretion (Bonnell, Vine 
& Bate 2004) and the link between the mass of a cluster and 
the mass of its most massive member (Weidner et al. 1201 Oi l. 
The distance is also of fundamental importance to understand 
the link with the TeV source HESS J1023-575 (Aharonian et 
al. 120071 1. which is spatially consistent with the bright Fermi- 
LAT pulsar PSR J1023.0-5746 (Abdo et al. [20l0l Ackermann 
et al. 1201 1| ). At much longer wavelengths, a peculiar jet as well 
as arc-like molecular ( I2 CO) emission feature were reported by 
Fukui et al. (j2009). The arc is well centred on the TeV source 
HESS J 1023-575, and Fukui et al. ( [20091 argue that both the 
TeV and GeV y-ray sources and the molecular features are likely 
to be associated with the same high-energy event. These au- 
thors accordingly suggest that these features could be due to 
an anisotropic supernova explosion of one of the most massive 
members of Westerlund 2 or to remnant objects such as a pulsar 
wind nebula or a microquasar. On the basis of the detection of 
significant emission during the off-pulse intervals with Fermi- 
LAT, Ackermann et al. ( 1201 11 ) argued in favour of the pulsar 
wind nebula scenario. In either case, the association of such a 
remnant with Westerlund 2 would imply that the most massive 
members of the cluster have already gone supernova, assum- 
ing star formation in this region was coeval. It must be stressed 
however that the preferred distance of PSR J1023. 0-5746 is 
2.4 kpc (Ackermann et al. 1201 11 ). Therefore, if the distance of 
Westerlund 2 were established to be much larger than 2.5 kpc, 
this would cast serious doubt on a connection between the pul- 
sar and the cluster. 

In this paper, we present the results of a spectroscopic 
follow-up campaign of the eclipsing binaries and the stars in 
the core of Westerlund 2. This campaign was specifically de- 
signed to address the questions outlined above. Throughout this 
work, we use the naming convention introduced by Moffat et al. 
( 11991b . except for objects that are discussed here for the first 
time. This paper is organized as follows. In Sect. [2] we provide 
a brief overview of our observing campaign and the data reduc- 
tion. Section[3]presents our results concerning spectral classifica- 
tion, variability, and multiplicity for each star. The implications 
of these results are finally discussed in Sect. [4] 



2. Observations 

We obtained eight observations of the newly discovered eclips- 
ing binaries in Westerlund 2 (Paper I), as well as several O- 
star binary candidates with the FLAMES multi-object spectro- 
graph on the Very Large Telescope (VLT) UT2 (see Table[T]). 
The observations were performed in service mode over an in- 
terval of one week in March 2008. Each observation lasted for 
an integration time of 2775 s and the observing conditions were 
good with a seeing of better than 1.1 arcsec. FLAMES offers the 
possibility to simultaneously acquire high-resolution spectra of 
up to eight targets with the UVES spectrograph and medium- 
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Fig. 1. Illustration of the position of the ARGUS unit (box) on 
the core of Westerlund 2. The top image uses our combined B 
and V image from Paper I (obtained with the SMARTS 1.3 m 
telescope at CTIO), whilst the bottom image uses an archive 
WFPC2 HST (F555W) image. The various objects analysed in 
the ARGUS data are labeled in the bottom figure. 



resolution spectra of some additional targets with the GIRAFFE 
spectrograph. 

Our campaign was divided into two parts (run A and B) 
with two different instrumental set-ups for the GIRAFFE spec- 
trograph. During all our observations, five well-isolated tar- 
gets (WR20b, MSP 18, as well as the three eclipsing systems 
MSP 44, 96, and 223) were observed with FLAMES-UVES fi- 
bres. The UVES echelle spectrograph was operated in the stan- 
dard setup Red 580, providing coverage of the wavelength do- 
mains 4842 - 5728 A and 5843 - 6800 A with a resolving 
power of 47000. During run A (five observations), the GIRAFFE 
spectrograph was fed by three integral field units (IFU), allow- 
ing us to monitor three additional O-type stars (MSP 171, 182, 
and 203). Whilst MSP 182 and 171 are well-isolated targets, 
MSP 203 lies in a more crowded region near the cluster core and 
its spectrum is heavily blended with that of MSP 444 (see Moffat 
et al. 119911 ). During run B (three observations), GIRAFFE was 
used with the ARGUS system. ARGUS consists of an array of 
14 x 22 microlenses offering integral field spectroscopy over a 
field of 11.8 x 7.3 arcsec with a sampling of 0.52". The field 
covered by ARGUS was oriented north-east and was centered 
at the position a =10:24:02.5, 6 = -57:45:24.0 (see Fig.[l} 
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Table 1. Overview of the targets studied in this paper. 



Object 


Inst. 


Epochs 


Binary status 


Previous SpT 


New SpT 


WR20b 


UVES 


8 


candidate 


WN6ha 


WN6ha 


MSP 18 


UVES 


8 


candidate 


05-5.5 V-m((f)) 


05-5.5V-ni((f)) 


MSP 44 


UVES 


8 


ecliping 


09.5 V 


Bl V + PMS? 


MSP 96 


UVES 


8 


eclipsing 


OB 


Bl: + Bl: 


MSP 167 


GIRAFFE-ARGUS 


3 


unknown 


06 III 


earlier than 06 


MSP 171 


GIRAFFE-IFU 


5 


unknown 


04-5 


04-5 


MSP 182 


GIRAFFE-IFU 


5 


unknown 


04 


04 


MSP 183 


GIRAFFE-ARGUS 


3 


unknown 


03V((f)) 


03V((f)) 


MSP 199 


GIRAFFE-ARGUS 


3 


unknown 


03-4 V 


03-4 V 


MSP 203+444 


GIRAFFE-IFU 


5 


unknown 


06v-m 


04-5 V 


MSP 223 


UVES 


8 


eclipsing? 


OB 


07-8 


A 


GIRAFFE-ARGUS 


3 


unknown 




08 V 


B 


GIRAFFE-ARGUS 


3 


unknown 




06-7 


C 


GIRAFFE-ARGUS 


3 


unknown 




07 V 


D 


GIRAFFE-ARGUS 


3 


unknown 




09.5: 


E 


GIRAFFE-ARGUS 


3 


unknown 




06-7: 



Notes. The first column lists the designation of the star, the second and third columns yield the instrument and the number of epochs for each 
target. The last three columns provide the binary status and spectral types prior to our study and the current spectral type (determined in this paper). 



Table 2. Radial velocities of prominent interstellar absorption and nebular emission lines in the FLAMES-UVES spectra of our 
targets. 



Line 


WR20b 


MSP 18 


MSP 44 


MSP 96 


MSP 223 




(kms- 1 ) 


(km s" 1 ) 


(kms" 1 ) 


(km s" 1 ) 


(kms" 1 ) 


Na i A 5890 


3.4 ±0.3 


2.8 ±0.3 


2.8 ± 0.5 


2.2 ±0.2 


2.3 ± 0.4 


Na i A 5896 


3.7 ± 0.2 


2.0 ±0.5 


2.9 ± 0.5 


1.8 ±0.3 


3.2 ±0.2 


[Om] ,14959 


17.1 ±7.6 


-3.7 ± 1.0 


28.3 ± 1.7 


22.3 ± 0.3 


7.4 ± 2.6 


[O m] ,15007 


18.6 ± 1.7 


-3.1 ± 1.5 


29.7 ± 1.2 


22.6 ±0.8 


8.2 ±0.9 



Notes. The error bars correspond to the 1-tr dispersion in our measurements 



near MSP 183 and MSP 16"Q For both runs (A and B), the 
GIRAFFE spectrograph was used in the low-resolution standard 
setup LR04 in third order covering the domain 5015 - 5815 A 
with a measured resolving power of 9500. 

All the data were reduced with the dedicated pipelines (ver- 
sions 4.2.1 for UVES, 2.5.1 for GIRAFFE-IFU and 2.8.1 for 
GIRAFFE-ARGUS). The signal-to-noise ratio (S/N) of the in- 
dividual FLAMES-UVES spectra reached between about 50 for 
the faintest stars and around 250 for WR20b and MSP 18. For 
the three FLAMES -GIRAFFE spectra, the S/N was around 150 
for MSP 171 and 182 and near 300 for MSP 203 + MSP 444. For 
the ARGUS data, the wavelength solution was obtained in two 
iterations yielding an rms in the calibration of 2.7 kms-'. The 
spectra were extracted using the SUM option. 

Owing to the high extinction towards Westerlund 2, the spec- 
tra of our targets are very rich in interstellar features. Besides the 
strong Na i AA 5890, 5896 absorption lines, we observe a large 
number of diffuse interstellar bands (DIBs, see Herbig 1995]). In 
addition, the spectra display a number of nebular emission lines 
from the RCW49 complex, such as [Om] AA4959, 5007 as weU 
as H i and He i. 

As can be seen from the results quoted in Table|2] the inter- 
stellar absorption lines are consistently found to indicate radial 
velocities (RVs) of around 2 - 3 kms" 1 for each star, whilst the 



2 The coordinates quoted in this paper are those from the archive 
WFPC2 HST image (based on GSC). They are shifted by -0.027 s in 
RA and -0.69" in Dec with respect to the 2MASS system used in Paper 
I. 



radial velocities of the nebular emission lines span a range from 
-3.7 to +29.7 kms" 1 , implying that the RCW49 nebula has a 
rather complex velocity field. 




532 534 536 538 540 



HJD - 2454000.0 

Fig. 2. Radial velocity variations about the mean evaluated from 
the five measurements for several sets of lines in the spectrum of 
WR20b (see text). The RV measurements of the narrow nebular 
[Om] ,15007 emission line are shown for comparison. 
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Fig. 3. Mean spectrum and temporal variance spectrum as computed from our time series of FLAMES-UVES spectra of WR20b. 
The dashed lines yield the 99% significance level of the TVS 1 '' 2 . 



3. Results 

3.1. WR20b 



WR 20b was reported as an Ha emission star by The ( I1966I his 
object THA 35-11-37). This result was confirmed by Schwartz, 
Persson & Hamann (119901 1. and the star was first classified as 
a WN7 Wolf-Rayet object by Shara et al. ( fT99TT >. This classi- 
fication was slightly revised to WN7:h, i.e. a WN7 star with a 
definite signature of hydrogen in its spectrum, by Shara et al. 
(1999), whilst van der Hucht (2001) proposed instead a spec- 
tral type WN6ha, i.e. a WN6 star with hydrogen emission lines 
and absorption lines belonging to the WN6 star. Subsequent ra- 
dio observations revealed the presence of a wind-blown shell 
around WR20b (Whiteoak & Uchida [T9971 ). Whilst the mean 
photometric properties of WR 20b do not support the presence 
of a luminous companion, the X-ray spectrum of the star is re- 
markably hard, a feature that is often associated with interacting 
wind binaries (Naze et al. 120081) . 

To address the issue of multiplicity, we measured the RVs of 
a number of stellar emission lines (H/?, Ha, He n AA 5412, 6406, 



6683, N iv A/15204, 5737, 6383, and N v A 4945, see Table|AT}, 
as well as interstellar and nebular features such as the diffuse in- 
terstell ar band s (DIBs) at AA5363, 5705, 5850, 6010, 6614 (see 
Herbig [19951 ), the Nai ^5890, 5896 doublet, and the [Om] 
A 5007 line. The interstellar and nebular features allow us to 
check the stability of the wavelength calibration: their 1 - <x RV 
dispersion is 0.3 km s for the narrow Na i absorptions, indicat- 
ing a very stable calibration. For the broad DIBs, the dispersion 
is larger, of the order of 4kms _1 . For the majority of the stellar 
emission lines, the RV dispersions are significantly larger than 
these values (24.8, 17.7, and 6.5 km s _I for the H i, He n and N iv 
+ v lines, respectively). However, as shown in Fig.[2] the varia- 
tions in the different sets of lines are uncorrelated. 

To quantify the spectral variability of the star, we computed 
the temporal variance spectrum (TVS, see Fullerton, Gies & 
Bolton [T9961 ) of our time series. This method provides a quan- 
titative assessment of the significance of the spectral line profile 
variability properly accounting for the noise level of the indi- 
vidual spectra. The results are shown in Fig. [3] Significant vari- 
ations are detected in the strongest emission lines, but not in 
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6500 6550 6600 

A (A) 



Fig. 4. Illustration of the Ha line profile in the spectrum of 
WR 20b for those two observations corresponding to the max- 
imum difference in apparent RV. The RV variations are clearly 
due to the top part of the line profile only. 



weaker ones such as the Niv AA 62 12-20 blend. The TVS of the 
strongest lines is quite complex, with multiple (up to four) peaks, 
indicating that the RV variability of these lines reflects their line 
profile variability. A more detailed inspection of the profiles of 
the hydrogen lines in different observations (see Fig. [4]) actu- 
ally indicates that the apparent RV variations are mainly due to 
changes in the upper part of the line profiles. The results pre- 
sented in this section suggest that WR 20b is indeed a spectro- 
scopically variable object, but this variability is more likely to be 
related to (large-scale) structures in the stellar wind, than caused 
by genuine RV variations of the star. Therefore, we conclude 
that there is at present no firm evidence that WR 20b is a short 
period (less than two week period) spectroscopic binary system 
(see also the discussion in Sect.[4|. 

Following the classification criteria of Smith, Shara & 
Moffat ( 1996), our spectra correspond most closely to a WN5- 
6 spectral type (ratio of the Hen A 5412 to Hei A 5876 intensi- 
ties above continuum = 3.6) with narrow lines (equivalent width, 
EW, of He ii A 5412 = -5.9 A > -40 A) and a definite signature 
of hydrogen emission revealed by the oscillation of the Pickering 
decrement (compare the intensities of Ha, Hen A 5412, and H/3 
on Fig[3]l. The only obvious stellar absorption feature in our 
spectra is the P-Cygni-type absorption trough of the He i A 5876 
line, which is obviously associated with the Wolf-Rayet star. 
Therefore, our data are consistent with the WN6ha classification 
given by van der Hucht (12001 I t. 

3.2. MSP 18 

On the basis of low -resolution spectroscopy, Moffat et al. d 1 99 1 I t 
classified MSP 18 as an 07:V star. This classification was re- 
vised to 04 V((f)) by Uzpen et al. (2005) and more recently to 
05-5.5 V-III((f)) based on higher resolution spectra (Paper I). 
Both its optical brightness (Paper I) and X-ray emission proper- 
ties (brightness, spectral hardness and slight flux variability; see 



-type stars in Westeiiund 2 
Table 3. Radial velocities of MSP 96. 



Date 


<Pi> 


4>RV 


Primary 


Secondary 


HJD-2454000 






(kms ') 


(kms ') 


532.587 


0.044 


0.058 




154 


532.631 


0.085 


0.099 


-182 


159 


534.581 


0.903 


0.917 


239 


-113 


534.628 


0.946 


0.960 


112 


-80 


536.624 


0.807 


0.821 


289 


-230 


537.626 


0.741 


0.755 


374 


-259 


538.719 


0.760 


0.774 


352 


-219 


539.622 


0.602 


0.616 


195 


-167 



Notes. The RVs are inferred from the multi-Gaussian deblending of 
the He i A 6678 line accounting for a nebular emission line. The second 
and third columns yield the orbital phases computed according to the 
photometric and spectroscopic ephemerides, respectively. 



Naze et al. 2008 ) suggest that this star could be an interacting- 
wind binary system. 

We fitted Gaussian line profiles to measure the RVs of a 
number of stellar absorption lines (H/3, Ha, Hei A 5876, Hen 
AA 5412, 6683), as well as two faint emission lines (C m AA 5696, 
6729-6731; Walborn|200D- The results are listed in Table[A2] 
We also measured the same interstellar and nebular features as 
for WR20b as well as the nebular H/3 and Ha emission lines. 
The 1 — cr RV dispersions of the stellar absorption lines are be- 
tween 1.8 and 5.9 kms -1 , which is quite comparable to the corre- 
sponding values of the prominent interstellar features. The stel- 
lar emission lines have slightly larger RV dispersions, but these 
emission lines are actually extremely weak (normalized intensity 
of 0.02 above the continuum for C in A 5696) and their RVs are 
thus subject to larger uncertainties. In addition, the TVS does not 
reveal any significant line profile variability. We thus conclude 
that MSP 18 is unlikely to be a short period (shorter than two 
week) spectroscopic binary (see also the discussion in Sect.Q. 

3.3. MSP 96 

The star MSP 96 was found to be an OB-type eclipsing binary 
system with a period of (1.0728 + 0.0006) days with probably at 
least one component close to filling its Roche lobe (Paper I). Our 
spectra do not reveal any trace of the He n A 5412 line, thus they 
are most consistent with a spectral type later than B0. 2 (Walborn 
1980). As a result of the severe contamination of the spectra by 
strong and asymmetric nebular emission lines in the hydrogen 
and neutral helium transitions, measuring the stellar components 
in the spectra of MSP 96 is quite a challenging task. 

The stellar spectral lines are broad, as expected for rapidly 
rotating stars in a contact binary system. Using a multi-Gaussian 
fitting technique, we attempted to deblend the Ha, He i AA 587 1 , 
6678 lines, accounting for the presence of a strong nebular emis- 
sion line. Our most reliable results were obtained for the Hei 
A 6678 line and the corresponding RVs are listed in Table[3] 
Unfortunately, the uncertainties on these RVs are quite large, of 
the order of 20 km s _1 . The He i lines are essentially identical in 
strength for the spectra of both stars, whilst the Ha line appears 
to be up to four times stronger in the spectrum of the star that has 
negative velocities during most of our observations. We caution 
however that this latter result is probably an overestimate be- 
cause of a bias resulting from the presence of the strong nebular 
emission line. 

In our data, the two components are most clearly deblended 
for the fifth, sixth, and seventh observations. Extrapolating the 



5 



Rauw, Sana & Naze: Early-type stars in Westeiiund 2 



X 



r 1.5 - 



— 

0) 
N 

a : 

o 

z 0.5 - 

0.04 F 

.03 h 

.02 =- 

.01 h 










DIB 








— 1 — 1 — 1 — 











> 

£ o 

^ 



4800 



5000 



5200 5400 
Wavelength (A) 



5600 




o 

y ' 0.5 - 

0.08 F 

> 0.06 =- 

E£ 0.04 : 

^ 0.02 E" I 









DIB 






H — 


— 1 — 1 — 1 — 















6000 6200 6400 6600 



Wavelength (A) 



Fig. 5. Same as Fig.^ but for MSP 18. 



light curve ephemerides of Paper I to the epoch of our spec- 
troscopic observations, we find that these observations should 
correspond to orbital phases near <p p = 0.75 (where <p p = 0.0 
corresponds to the primary eclipse, i.e. the brighter star being 
eclipsed). Although our spectroscopic data are taken 1162 days 
after Tq of Paper I and despite the short orbital period, the or- 
bital solution below agrees with the photometric ephemerides. 
Indeed, the shift in phase between the photometric and spectro- 
scopic solution is found to be 0.014. 

We note that the star with the stronger Ho- line is seen 
blueshifted at phases near <p p = 0.75, hence that this line be- 
longs to the intrinsically fainter component (i.e. the photometric 
secondary) of the system. From these considerations, we con- 
clude that the photometric secondary must be intrinsically some- 
what cooler than the primary, but should have a somewhat larger 
radius. The helium line intensity ratio might then be close to 
unity, although the primary star (i.e. the one that is eclipsed at 
phase 0.0) would be somewhat hotter, thus explaining the deeper 
primary eclipse. Alternatively, the star with the negative RVs 
near <p p = 0.75 could be the primary star, but this would im- 
ply a shift in phase between the photometric and spectroscopic 



ephemerides of 0.486. This seems a priori rather unlikely, but 
cannot entirely be ruled out given the error in the orbital period 
and the large number of orbital cycles between our photometric 
and spectroscopic campaigns. We briefly consider the implica- 
tions of this alternative assumption below. 

3.3.1. Orbital solution 

To establish the orbital solution of MSP 96 from the RV data 
in Table[3] we started by considering the SB1 solution for the 
photometric secondary star. For this purpose, we used the Liege 
Orbital Solution Package (LOSP) code (Sana & Gosset l2TjTTI 
which allowed us to derive K s = (270 + 22)kms~ 1 , as well as 
y = (19 + 19)kms _1 . We then considered the RVs of the photo- 
metric primary star with the value of y fixed to 19 km s _1 , which 
yielded K p = (337 ± 78)kms~ 1 . From these results, we then in- 
ferred a secondary/primary mass ratio of 1 .25 + 0.27, as well as 
m p sin 3 / = (11.1 ± 3.0) M G and m s sin 3 / = (13.8 ± 6.6) M B . 
These minimum masses, along with the eclipsing nature of the 
system, which indicates a rather large inclination, suggest a spec- 
tral type near Bl for both components. This classification is also 
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consistent with the overall spectral morphology, although the 
large number of interstellar features and the moderate S/N of 
our data prevent us from being more specific. 



3.3.2. Light curve solution 

Adopting the mass ratio inferred hereabove, we then analysed 
the photometric data of the system from Paper I using the night- 
FALirl code, developed and maintained by R. Wichmann, M. 
Kuster and P. Risse. We set the temperature of the primary star 
to 25400 K (appropriate for a B 1 V classification, Schmidt-Kaler 
1982 ) and allowed the inclination, the Roche-lobe filling factors 
(hereafter fill p and fill s ), and the secondary temperature to vary. 
The best-fit model parameters were found to be i = (76.8^ 2)°, 
T s = 21550+™ K, fill,, = 0.84 ± 0.03, mid fill, = 0.85 ± 0.02. 
This inclination yields absolute masses and radii, respectively 



of M p = (12.0 ± 3.3) Mq, M s = (15.0 ± 7.2) M Q , and R p = 
(3.8 ± 0.6)R Q , R s = (4.3 + 0.6) Rq. These radii are rather small 
for the masses, but we note that quite similar parameters were 
found for other early B-type main-sequence stars in eclipsing 
binaries (e.g. V 578 Mon and DW Car, see Southworth |2007l and 
references therein). That the more massive star (the secondary) 
would be the cooler and less luminous component of the sys- 
tem is surprising and indicates that this star would actually be 
more evolved. The radii of both stars suggest instead that they 
are rather unevolved main-sequence stars. Whilst a case A bi- 
nary mass transfer could alter the properties of the components, 
we stress that the light curve solution relates to stars that do not 
fill up their Roche lobes. This issue would actually be resolved 
if we assumed a phase shift of 0.486 between the photometric 
and spectroscopic ephemerides (see below). 




Fig. 6. Combined solution of the light curve and the RV curve of 
the MSP 96 system (see text). The orbital phases are computed 
from the photometric ephemerides of Paper I. 



Combining all the above parameters, we inferred a total 
bolometric luminosity of the binary of 8960 L G . The uncertainty 
in this number is difficult to evaluate, because of the unknown 
uncertainty in the primary's spectral type. However, we esti- 
mated that the relative uncertainty is likely to be about 25%. 
From V = 16.17 at phases outside the eclipses and E(B — V) = 
1 .50 as follows from Fig. 7 of Paper I, we then inferred a distance 
of 6.6 kpc again probably subject to a 25% uncertainty. 

We also analysed the light curve assuming a shift of 0.486 in 
phase between the photometric and spectroscopic ephemerides. 
In this case, the mass ratio inferred in Sect. |3.3.l"| would be in- 
verted, m s /m p = 0.80 + 0.17, and the more massive star (which 
is now the primary with m p sin/ = 13.8 M ) would also be 
the brighter one. The corresponding best-fit photometric solu- 
tion is essentially symmetric with respect to the parameters in- 
dicated above (z = 77.2°, T p = 25400 K (fixed), T s = 21790 K, 
fillp = 0.85, fill s = 0.84) and yields a similar distance of 6.8 kpc. 
This solution corresponds to an unevolved B + B main-sequence 
binary system and avoids the problems described above. 

3.4. MSP 44 

In Paper I, we classified MSP 44 as an eclipsing binary with a pe- 
riod of 5. 176 days, noting that the aliasing period at 1 .2395 days 
yielded an almost equally strong peak in the periodogram. 
However, the rather narrow eclipses argue against the shorter or- 
bital period and we thus favour the solution corresponding to 
the (5.176 + 0.029) day period. Extrapolating the photometric 
ephemerides from Paper I to the epoch of our FLAMES observa- 
tions, we find the phases listed in the second column of Table|4] 
where phase 0.0 corresponds to the primary eclipse. 

Whilst the primary eclipse has a depth of 0.6 mag, the sec- 
ondary eclipse was found to be much shallower (< 0.1 mag). 
There are two possible interpretations of this situation: either 
the surface brightnesses of the two components of the system 
are very different or the binary orbit is sufficiently eccentric with 
a longitude of periastron (for the star being eclipsed at phase 0.0) 
near 90f\ In addition, the light curve exhibits some slow modu- 
lation, especially around secondary minimum that we tentatively 
attributed to spots on one of the stars (Paper I). 

The absence of any Hen line in the spectrum of MSP 44 
clearly indicates a spectral type later than B0. 2 (Walborn 1980) 
for the combined spectrum of MSP 44. This is slightly at odds 
with the 09.5V: photometric spectral classification of this star 
(Naze et al. 2008 ). The likely presence of a very weak Si in 
A 5740 feature and possibly of Si n A 6347 suggest a spectral type 
near Bl with an uncertainty of at least one subtype. Given the 
probably rather large brightness ratio of the primary to the sec- 
ondary, we can actually adopt this spectral type as the one of the 
primary component. The Bl spectral classification is also sur- 
prising in view of the rather strong and hard X-ray emission of 
MSP 44 measured with Chandra (Naze et al. l20081 >. 

Measuring the stellar RVs was, once more, quite difficult. 
The reasons for this were (1) the contamination of most, suffi- 
ciently intense, stellar lines (H i and He i) by rather strong neb- 
ular emission lines, and (2) the moderate S/N of our data that 



4 If the eccentricity were sufficiently large and the longitude of pe- 
riastron close to 90°, one could expect to observe a difference in the 
durations of the primary and secondary eclipses, although the ampli- 
tude of this effect strongly depends on a number of parameters (orbital 
eccentricity and inclination, filling factors,...) and is sometimes even 
absent. Within the accuracy of our data, no such effect is present in the 



http://www.hs.uni-hamburg.de/DE/Ins/per/Wichmann/Nightfall.html light curve of MSP 44. 
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Table 4. Radial velocities of MSP 44. 
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536.624 
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0.247 
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29.6 


539.622 


0.632 


0.516 


45.1 


-264 


315 


27.6 



Notes. The second and third columns yield the orbital phases computed according to the photometric and spectroscopic ephemerides respectively. 
RVi , RV A , and RV r stand for the RVs of the main component of the He i A 5876 line and the blue and red components of this line, respectively (see 
text). 



hampers the detection of very faint spectral lines. At first sight, 
only the signature of one star is seen in the spectra. The RVs 
of this main component vary very little during our observations 
(which cover a priori both quadrature phases of the orbit). This 
is very puzzling and prompted us to re-inspect the morphology 
of several lines in more detail. 

The most intriguing line turned out to be Hei /I5876. The 
blue and red wings of this line suggest the presence of two addi- 
tional absorption features in most of our data (see Fig.|7}. There 
are no known DIBs at these wavelengths (Herbig 1995 Tjand the 
blue wing of the line at least should not be affected by any strong 
telluric absorption lines either (Curcio, Drummeter & Knestrick 
1964). 
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Fig. 7. The region around the He i A 5876 line in the spectrum of 
MSP 44 observed on HJD 2454539.622. The emission peak near 
5876 A is the nebular He i line. Note the complex morphology of 
the line, especially in the blue and red wing. 



Since it is possible that MSP 44 might be a triple system, 
we thus performed a multiple Gaussian fit of this line with up 
to four components (three stellar absorption components and 
one nebular emission component). The nebular component was 
consistently fitted with an RV of (29.4 + 3.1)kms _1 . The main 
stellar component was in turn found to display RVs between 
17 and 46 kms 1 . The corresponding RVs are listed in Tableji] 
Although the RVs of the blue- and red-shifted absorption fea- 
tures must be interpreted with extreme caution, it is immediately 



An up-to-date catalog of DIBs can be found 



clear that if MSP 44 were a triple system, then the mass ratio of 
the binary responsible for the blue and red absorption compo- 
nents would be close to unity. This in turn would be at odds with 
the light curve, which implies that there is a rather large differ- 
ence in surface brightness between the components of the eclips- 
ing binary. The most likely interpretation of this line morphol- 
ogy therefore seems that it represents the very broad absorption 
line of only one star, but either deformed by pulsations such as in 
/3Cep-type variables (e.g. Telting et al. 2006) or partially filled in 
by some circumstellar emission or chromospheric emission from 
a low-mass companion heated by the irradiation from the B-star. 
In the case of the former interpretation, it must be stressed that 
our photometry does not reveal the signature of typical ySCep 
variations (see Paper I). As to the circumstellar emission inter- 
pretation, we note that neither H/3 nor Ha have emission above 
the continuum, but both lines have a somewhat similar morphol- 
ogy to that of He i A 5876, in the sense that they exhibit a rather 
narrow core and more extended wings (especially the red wing). 

In the absence of an unambiguous interpretation of the spec- 
tral morphology of MSP 44, we choose to use the RVs of the He i 
A 4921 line, which is less complex than Hei /15876, to derive 
some constraints on the orbital motion of the main component. 
Assuming that the RVs of this line indeed reflect orbital motion, 
a circular SB1 orbital solution yields T = 2454531.774 + 0.285, 
71 = (42+9) kms" 1 , K x = (34+ 16) kms" 1 , and a x sin/ = 3.5tf . 

This yields a very low mass-function of f(m) 



at 



http : //leonid . arc . nasa . gov/DIBcatalog . html 



(mi+m2) 

(0.021+0.030) M . Here, To stands for the time of primary mini- 
mum. Compared to the photometric ephemerides, there is a shift 
in phase by 0.116 or 1.116. Since the number of orbital cycles 
between our photometric and spectroscopic campaigns amounts 
to 224, this phase shift implies an error in the orbital period of 
either 0.003 or 0.026 days, well within or consistent with the es- 
timated error in the orbital period as inferred from the photomet- 
ric data (0.029 days, see above). 

Since the system is eclipsing, we can to first order approxi- 
mate the sin 3 i term in the mass function to be unity. Assuming 
a mass of 10 - 12 M Q for the Bl primary, we find that the sec- 
ondary should have a mass of about 1.4 - 1.6 M G , which implies 
that the companion might be an early F-type star. But would this 
result be reasonable? To answer this question, we return to the 
photometric light curve of the system. We first recall that the lat- 
ter displays variations outside the eclipses that are quite large. In 
Paper I, we attributed them to spots, but our current knowledge 
of the system allows us to be a bit more specific. If we are indeed 
dealing with an F star closely orbiting around a B 1 V primary, it 
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seems unavoidable that there will be a strong heating effect of 
part of the secondary's surface by the primary's radiation. Test 
calculations carried out with the nightfall code indicate that the 
light curve of the system can a priori be explained by assum- 
ing that the hemisphere of the secondary facing the primary is 
significantly hotter than the rear side. However, the solution of 
the light curve is certainly not unique, because we ignore the ac- 
tual mass ratio of the system and this parameter cannot be well- 
constrained from the photometric data alone. As an illustration, 
we show below the best-fit solution obtained for mi/mj = 8. 




> 



Fig. 8. Combined solution of the light curve and the SB1 RV 
curve of MSP 44. The parameters of the light curve model are 
1112/1111 = 0.125, i = 80.6°, primary and secondary radii of 
5.2 and 5.6 R G respectively, and mean surface temperatures of 
Ti = 29100 K (fixed) and T 2 = 6440 K. The secondary star fea- 
tures a hot spot facing the primary with a radius of 45° and a 
temperature enhanced by a factor 1 .76. 



We note in addition that the radius of the secondary inferred 
from the light curve solution (5.6 R Q in the case of the best-fit 
solution shown below) is actually much larger than the expected 
radius of an early-F main-sequence star, but closer to what one 
would expect for a giant. Given the youth of the Westerlund2 
cluster, the most straightforward explanation is that the sec- 
ondary star of MSP 44 is not a genuine giant, but rather a pre- 
main sequence (PMS) object. A PMS nature of the secondary 
would also provide a very natural explanation of the rather bright 
and hard X-ray emission of this system. 

In summary, our analysis of the spectroscopic and photomet- 
ric data of MSP 44 leaves a number of open questions about the 
exact nature of the primary component. The most likely scenario 
seems to be that this system consists of a B 1 main-sequence star 
and a low-mass PMS companion. The physical parameters of the 
system are admittedly too uncertain to attempt a self-consistent 
determination of the distance of Westerlund 2, although we note 
that the parameters of the fit in Fig. [8] are consistent with a dis- 
tance of 8.4 kpc. If the Bl V spectral type of the MSP 44 pri- 
mary were correct and the companion contributed little to the 



integrated light of the system, then the V magnitude and B — V 
colour outside eclipse would indeed be consistent with a distance 
near 8.4 kpc. 

3.5. MSP 223 

The yellow -red spectra of MSP 223 display a rather strong He 11 
A 5412 feature, as well as a weak (EW ~ 0.4 A), but definitely 
detected O in A 5592 line (see Fig.|9]l. Whilst there are no known 
interstellar features affecting the O m line, there are three, rather 
weak, DIBs at AA 5404.5, 5414.8, and 5420.2 that are blended 
with Hen A 5412. To assess the real strength of the Hen A 5412 
feature, we first evaluated the ratio of the EWs for some nearby 
DIBs (/U5450.3, 5487.5) in the spectrum of MSP 223 and the 
corresponding EWs in the spectrum of HD 183143 quoted by 
Herbig ( 119951 ). This ratio (1.33) was then used as a linear scaling 
factor to estimate the contamination of the Hen A 54 12 line by 
the three DIBs. The corrected EW of the He 11 line amounts to 
0.9 A. Therefore, Hen A 5412 is of comparable strength to the 
He 1 A 5876 line (EW = 0.85 A), although the latter is slightly 
affected (probably less than 0.2 A) by a nebular emission line. 
These properties suggest an 07-8 spectral type (Walborn 1980) 
for the combined spectrum. That we are dealing with an O-type 
star is confirmed by the absence of Hei A 4921 and the fact that 
the He 1 A 6678 line is clearly blended with He 11 A 6683. 
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Fig. 9. Mean spectrum of MSP 223. Most of the features are due 
to interstellar lines or bands. The rather prominent Hen A 5412 
line indicates an O-star classification (see text for details). 



We measured the RVs of several stellar lines. The most ro- 
bust results were obtained for Hen A 5412 and Ha, although we 
caution that the former is affected by blends with DIBs and for 
the latter, we had to perform a two Gaussian fit to account for 
the contamination by the nebular Ha emission line^] The results 
are listed in Tableland plotted in Fig. 10 after subtracting the 
mean RV individually for each line. No obvious trend is appar- 
ent in this figure. The rather large shift in RV between the He 11 



In these fits, the nebular line has a very stable RV of 1 . 1 ±0.2 km s 
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Table 5. Radial velocities of MSP 223. 
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64.2 


-54.2 
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60.0 


-36.0 


539.622 


86.3 


-24.8 



A 5412 and Ha line is most likely due to wind emission affecting 
the latter line. Similar shifts in the apparent systemic velocity of 
different lines are commonly found in binary systems where at 
least one component features a strong stellar wind (e.g. Rauw et 
al. 120021 . 
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Fig. 10. Radial velocity variations of MSP 223 during our ob- 
servations. The error bars in the He RVs yield the formal error 
in the position of the absorption line in the two-Gaussian fit ac- 
counting for the nebular emission line. 

Our previous photometric monitoring of the cluster indi- 
cated that MSP 223 is probably an eclipsing binary (Paper I). 
The Fourier analysis yielded the highest peak for a period of 
14.61 days, with a rather large uncertainty of 1.2 days. The light 
curve folded with this period yields only a primary eclipse, with 
no indication of a secondary eclipse, although we might have 
missed it, if its duration is sufficiently short to fit into a small 
gap in phase coverage around (p = 0.5. In Paper I, we inferred 
an orbital period of 29.21 1 days, which then indicates two rather 
similar eclipses. This would imply that the system consists of 
two rather similar stars, well within their Roche lobe and seen at 
a large inclination angle. However, we do not observe an obvious 
SB2 signature in our spectra of MSP 223, which suggests instead 
that the companion of the O-star is significantly fainter and less 
massive. This would be consistent with the rather modest RV 
variations that we have measured (see Fig. [TO}, which show a 
peak-to-peak amplitude of about 40kms~ 1 over seven days (the 
most reliable RVs are probably those of Ha). Unfortunately, the 
considerable uncertainties in the orbital period prevent us from 
extrapolating the photometric ephemerides to the epoch of our 
spectroscopic observations and we cannot derive constraints on 
the orbital solution of this system. 

Although our knowledge about MSP 223 is still rather poor, 
we can nevertheless attempt to derive its spectro-photometric 
distance. For this purpose, we combine the photometric informa- 



tion from Paper I (V = 15.75, B -V — 1.50) with the absolute 
magnitudes and intrinsic colours of 07-8 V stars from Martins & 
Plez (2006 ). In this way, we obtain a distance modulus of 14.67 - 
14.96 corresponding to a distance of 8.6 - 9.8 kpc. If our spectral 
classification were too early, and we adopted instead an 09.5 V 
type, we would still derive a distance of 7.2 kpc. We stress that 
this number does not account for the contribution of a compan- 
ion to the light of MSP 223. Accounting for this would lead us 
to inferring an even larger distance. 

3.6. MSP 1 71, 182, and 203 + 444 

These three (actually four) objects were observed with the inte- 
gral field spectroscopy units connected to the GIRAFFE spec- 
trograph. Their mean spectra are shown in Fig.[TT| 
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Fig. 11. Mean GIRAFFE spectra of MSP 171, 182, and 203 + 
444. 



In Paper I, the combined blue spectrum of MSP 203 + 
MSP 444 was assigned a spectral type 06V-III. However, the 
GIRAFFE spectra of MSP 203 + MSP 444 did not reveal any 
trace of Hei AA5016, 5048. Therefore, by comparison with 
the yellow-green spectra of OB stars published by Kerton, 
Ballantyne & Martin (11999b . the brightest star of the blend must 
have a spectral type earlier than 06 V and possibly even earlier 
than 05 V. The mean spectrum of MSP 203 + MSP 444 is very 
similar to that of MSP 171 and 182, which were classified as 
04-5 and 04 respectively in Paper I. We thus conclude that the 
brighter component of MSP 203 + MSP44^]is likely of spectral 
type 04-5 V. The other component is likely responsible for the 
weak He i A 4471 absorption that was seen in Paper I. We note 
that the latter absorption is significantly broader than the other 
stellar absorption lines in the blue (and yellow-green) spectra of 
MSP 203 + MSP 444. This supports our suggestion that this line 
is not associated with the same star as the stronger, higher ion- 
ization lines. 

We measured the RVs of the three stars by fitting Gaussian 
profiles to the Hen A 5412, Om ^5592, and Civ /LI 5801, 5812 
lines with the caveat that the Civ lines (especially Civ A 58 12) 
are heavily blended with DIBs. We also measured the RV disper- 



7 According to the photometry of Moffat et al. J199U . MPS 203 is 
slightly brighter than MSP 444. However, the HST data (see Fig.[T]( sug- 
gest the reverse situation. 
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sions of four DIBs (AA 5362, 5487, 5705, and 5781, see Herbig 
119951 , which are of comparable strength to or slightly stronger 
than the stellar lines. 

For MSP 171, we infer a 1 - cr RV dispersion of 6.9 km s 1 
for the strongest stellar line (He n A 5412) and 3.5 km s _I for the 
Civ A 58 12 line. For the same set of spectra, the RV dispersions 
of the DIBs are between 0.8 and 10.3 km s _1 with a typical value 
of about 5kms _1 . 

In the case of MSP 182, all stellar lines, except for Civ 
,15812, which is affected by the DIB at 5809 A (Herbig [T995I 
display a 1 - cr RV dispersion in the range 3.0 - 5.7 km s -1 . For 
the interstellar bands, these dispersions are 1.1 - 4.1 kms _I , ex- 
cept for the rather weak DIB at A 5362 which has a larger disper- 
sion of 9.3kms~ 1 . 

Finally, the strongest stellar lines of the MSP 203 + MSP 444 
spectra have RV dispersions of 2.5 - 4.1 km s _1 , whilst the weak 
O in A 5592 feature has a slightly larger dispersion of 7.6 km s _1 . 
In this case, the RV dispersions of all DIBs measured amount to 
1.6-3.3 kms- 1 . 

In conclusion, none of the stars observed with the IFUs feed- 
ing GIRAFFE is found to display significant RV variations on 
timescales of a week. 

3.7. Stars observed with the ARGUS unit 

The field of view covered by the ARGUS unit harbours two stars 
with an MSP identifier (MSP 167 and 183) as well as five anony- 
mous objects (labeled A to E in Fig.[T]l that are bright enough to 
extract at least an averaged spectrum. In addition, one spaxel at 
the edge of the field contains photons from the wing of the PSF 
of star MSP 199 which is right outside the nominal field of view. 
For the brightest and most isolated objects, the spectra were ex- 
tracted over nine spaxels (MSP 167, MSP 183, object A), whilst 
the extraction region was reduced to a single spaxel for those 
stars (objects B, C, D, and E) that are either relatively close to a 
bright source or fall only partially inside the field of view of the 
ARGUS unit (MSP 199). In this section we discuss the results 
obtained for the different stars in the ARGUS field. 

In Paper I, MSP 183 was classified as an 03 V((f)) star. The 
lack of He i ,1,15015, 5047 absorption, the strength of the Hen 
A 5412 line and the presence of the N iv AA 5200 - 5205 blend in 



the GIRAFFE- ARGUS spectrum of MSP 183 (see Fig.[T2]i are 
in line with this classification (Kerton et al. 1999). However, the 
most remarkable feature of the spectrum of MSP 183 is the pres- 
ence of the Civ ,1,15801, 5812 lines in weak, but definite, emis- 
sion with equivalent widths of -0.24 A and -0.06 A, respec- 
tively. These lines were reported to be in emission in the spec- 
trum of the 03 If* supergiant Cyg OB2 #7 (Walborn & Howarth 
2000). Our blue spectrum of MSP 183 (see Paper I) rules out the 
possibility that MSP 183 could be such an evolved supergiant. 
Therefore, it seems more likely that these C iv emission lines are 
not selective properties of 03 If* supergiants, but could actually 
be a general feature of the hottest O-type stars as suggested by 
Walborn (1200Tb . 

The RVs of MSP 183 were measured by fitting Gaussian pro- 
files to the strongest stellar lines (Hen A 5412, Om ,15592, Civ 
AA 5801, 5812; see Table|A6]l. In this way, we obtain a 1 - cr RV 
dispersion of 1.2kms~' for the Hen A 5412 absorption and 7.3 
and 3.3 kms -1 for the Civ A 5801 and A 58 12 emission lines re- 
spectively. The RV dispersion of the DIBs range between 0.9 and 
7.5 kms 1 depending on the strength and shape of the DIB. We 
therefore conclude that this star does not display any significant 
RV variations during our observations. 
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Fig. 12. Mean GIRAFFE- ARGUS spectra of some of the targets 
discussed in Sect.13.71 



On the basis of a spectrum in the blue domain, we previously 
classified MSP 167 as an 06 III star (Paper I). The yellow spec- 
tral domain, however, is more consistent with an earlier spectral 
type. Indeed, our FLAMES-ARGUS spectra reveal no trace of 
either the He i AA 5015, 5047 or the O m A 5592 absorption lines. 
Furthermore, the C iv AA 5801, 5812 lines are clearly not present 
in absorption. If anything, Civ ,15801 is in very weak emission 
(EW 0.10 A). These properties are at odds with an 06 spec- 
tral type because the Om /15592 and C iv /LI 5801, 5812 absorp- 
tion lines are expected to be strongest for the spectral types 06- 
07 (Walborn[l9E)]). 

The dominant stellar line is Hen A 5412. We measured its 
RVs by fitting Gaussian profiles (see Table|6]l. Whilst there is 
only little difference between the RVs of the first two spectra, 
there is a clear shift of more than lOOkms -1 between the sec- 



ond and third observation (see Fig. 13 1. This shift is highly sig- 
nificant since the RV dispersion of the DIBs range between 1 .0 
and 9.5 km s with the lower limit on the dispersion holding for 
those DIBs which are of similar strength to the He n A 5412 line. 
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Fig. 13. The He n A 5412 line in the spectrum of MSP 167 at the 
three dates of our GIRAFFE- ARGUS observations. The shift be- 
tween the second and the third observation is clearly seen. 



The above results suggest that MSP 167 is actually a short- 
period (of the order of a week) spectroscopic binary consisting 
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Table 6. Radial velocities of the binary candidate MSP 167. 



HJD-2454000 


RV (He a A 54 12) 




(tans- 1 ) 


532.631 


-56.2 


534.628 


-64.5 


537.626 


49.0 



of an early O-type star (responsible for the strong He n line seen 
in our ARGUS data) and a companion of later spectral type 
(which produces the He i absorption lines seen in the blue-band 
spectra). We note that this star displays a rather large X-ray 
luminosity of 8.9 x 10 32 ergs~' with some hints of moderate 
(2 cr) variability (see Naze et al. 2008 ). These properties, along 
with the RV variations reported above, suggest that MSP 167 
might be an interacting-wind binary system. 

For MSP 199 (03-4 V according to Paper I), the spectra were 
extracted from a single spaxel of the ARGUS unit. Therefore, 
the S/N of the individual spectra is rather low. Nevertheless, the 
combined spectrum agrees with the classification proposed in 
Paper I. The RVs of the Hen A 5412 display a 1 - cr dispersion 
of 5.3 kms 1 , whilst the DIBs at ,15781 and ,15796 have RV 
dispersions of 1.9 and 2.4kms~ 1 , respectively. 

Objects A, B, C, D, and E are discussed here for the first 
time. Our mean spectrum of object A (10:24:02.8 -57:45:30.7, 
J2000) displays rather strong Hen ,15412, Om ,15592, and Civ 
,1,15801, 5812 absorption lines along with fainter absorption 
lines of Hei ,1,15016, 5047 and a very weak Cm ,15696 
emission line (see Fig.[T2"|). These spectral features are indicative 
of an 08 spectral classification with an uncertainty of one 
subtype (Kerton et al. 1999, Walborn 1980). The weakness 
of the Cm emission most likely indicates a main-sequence 
luminosity class. Our photometric data (Paper I) indicate 
V = 15.766 ± 0.050 and B - V - 1.456 for this star. Assuming 
Rv = 3.1 and adopting the intrinsic colours and absolute 
magnitudes of an 08 V star (with an uncertainty of one spectral 
type) from Martins & Plez (2006), we infer a distance modulus 
14.83 + 0.30 corresponding to a spectrophotometric distance of 
(9.2 + 1.3) kpc. The RVs of the various stellar absorption lines 
yield 1 - cr dispersions in the range 1.1 to6.5kms 1 depending 
on the strength of the line considered (see Table |A.8| . The 
lowest dispersions are found for the He n A 5412 and Civ A 5801 
lines. The DIBs at ,15781 and ,15796 have RV dispersions of 
4.1 and 1.3 km s , respectively. 

For star B (10:24:02.6 -57:45:31.2, J2000), the lines of Hei 
and O in are weak, but clearly present, whilst the C iv absorption 
lines are also rather weak (see Fig. 12 1. The strength of the He n 
,15412 line implies a spectral type around 06-7. The RVs of the 
latter line indicate a 1 - cr dispersion of 3.8kms _1 , which is 
slightly larger than the one of the DIB at ,15781 (1.9 kms -1 ). 

Star C (10:24:02.5 -57:45:32.1, J2000) displays a weak 
Hei ,15016 line, as well as shallow and broad Om and Civ 
absorption lines (Fig. 12 1, and Cm ,15696 is seen in very 



weak emission. The spectral type is likely about 07 V. The RV 
dispersion of the Hen ,15412 line amounts to 2.2kms~ 1 , very 
much comparable to that of the DIBs at A 578 1 and A 5797 (2.3 
and 1.6kms~', respectively). 



For stars D and E, the S/N of the spectra is rather low, 
preventing us from measuring the RVs of the stars from the 
individual spectra. We instead combined the three observations 
of each star into a mean spectrum with the goal of performing 
a spectral classification. Object D (10:24:02.2 -57:45:32.0, 
J2000) has a stronger Hei ,15016 than Hen ,154 12 line . We 
accordingly classify this star as 09.5 (Kerton et al. 1999). Star 
E (10:24:02.6 -57:45:32.9, J2000) must be of earlier spectral 
type, since there is no trace of Hei ,15016, whilst Hen ,15412 
is the strongest stellar line in the spectrum. A weak O m A 5592 
line is seen, but no Civ absorption could be detected. We 
suggest a spectral type near 06-7. 

We inspected the Chandra image of the cluster's core (see 
Naze et al. 120081 1 at the positions of objects A - E. Star E cor- 
responds to the southwest lobe of an X-ray source detected be- 
tween MSP 183 and MSP 199 using the wavedetect algorithm 
(Naze et al. 2008). Whilst there is no obvious X-ray emission 
from stars A and C, there might be some X-ray emission asso- 
ciated with stars B and D. However, the latter two objects lie 
in the wings of the X-ray source associated with MSP 183, and 
their status as X-ray emitters can therefore not be ascertained 
with 100% confidence. The lack of X-ray detections for stars A 
and C is somewhat surprising given their spectral types inferred 
above. One possibility could be a larger interstellar column den- 
sity, thus larger absorption of the X-ray emission than for other 
stars in Westerlund2 of similar spectral type. However, this ex- 
planation needs to be confirmed using additional high angular 
resolution optical and X-ray observations. 

4. Discussion 

The various eclipsing binaries that we have studied in this paper 
yield distance estimates in the range 6.5 to larger than 9 kpc. 
Whilst the uncertainty in these results is difficult to evaluate, 
they are in very good agreement with the spectrophotometric 
distance of (8.0 ± 1.4) kpc inferred in Paper I and clearly 
consistent with a cluster distance well beyond the 2.8 kpc 
proposed by Ascenso et al. (2007). We thus conclude that all 
available observations of the early-type stellar population of 
Westerlund2 are most consistent with a distance near 8.0kpc 
in agreement with the distance inferred for WR 20a. This result 
casts serious doubts on a possible association between the 
y-ray pulsar PSR J 1023. 0-5746 (Ackermann et al. [20TTb and 
the Westerlund2 cluster. If the distance of the pulsar were 
confirmed at 2.4 kpc, it would then appear more likely that the 
latter belongs to a foreground stellar population unrelated to the 
massive stars of Westerlund 2, but possibly related to the PMS 
stars discussed by Ascenso et al. (120071 1 . 

With respect to the multiplicity of early-type stars in 
Westerlund 2, we stress that among the eleven O-type stars and 
one WNha star monitored during our campaign, which were not 
previously known to be binary systems, only one star (MSP 167) 
was found to probably be a short-period spectroscopic binary 
system. The other targets did not display RV variations well 
above the uncertainties estimated from the RV dispersion of the 
DIBs. To assess the significance of our results, we performed 
Monte Carlo simulations (Hammersley & Handscomb 1 1 9641 l of 
the RVs of a synthetic population of massive binaries. The parent 
distribution of the binary parameters of the population of mas- 
sive binaries is unfortunately unknown, but we can make a series 
of reasonable assumptions. Here, we have followed an approach 
similar to the one of Sana, Gosset & Evans (120091) . Following 
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Table 7. Results of our Monte Carlo simulations of the detection of massive binaries with the sampling of the GIRAFFE-IFU 
campaign. 



mi (M B ) 


I 


II 


III 


IV 


V 




P< 5 


5 <P< 10 


10 < P < 20 


20 < P < 50 


50 < P < 100 




(%) 


(%) 


(%) 


(%) 


(%) 


20 


0.7 


0.7 


1.8 


19.0 


56.2 


30 


0.6 


0.6 


1.5 


15.8 


50.8 


40 


0.5 


0.5 


1.2 


13.8 


46.6 


50 


0.4 


0.4 


1.0 


12.4 


43.6 



Notes. For a given range of orbital periods (expressed in days) and a given value of the primary mass, the table yields the fraction of systems (in 
% of the population of systems within this interval of periods) that would display ARV = max RV(<pj) - min RV((f>j) < 20kms~', hence escape 
detection. The total numbers of simulated systems in each of the period intervals are 34198, 25802, 4786, 6315, and 4772 for the intervals I, II, 
III, IV, and V, respectively. 

Table 8. Same as Table^ but for the sampling of the GIRAFFE- ARGUS data set. 



m, (M ) 


I 


II 


III 


IV 


V 




P< 5 


5 < P < 10 


10 < P < 20 


20 < P < 50 


50 < P < 100 




(%) 


(%) 


(%) 


(%) 


(%) 


20 


3.4 


1.5 


4.5 


29.8 


70.8 


30 


2.7 


1.2 


3.4 


26.2 


64.8 


40 


2.4 


1.0 


2.8 


23.7 


60.7 


50 


2.1 


0.9 


2.3 


21.7 


57.6 



the latter authors, we adopt a bi-uniform period distribution in 
log P(days) with 60% of the systems in the range 0.3 < log P < 
1 .0 and the remainder in the range 1 .0 < log P < 3.5. The eccen- 
tricity was assumed to be uniformly distributed between 0.0 and 
0.9, with all systems with orbital periods shorter than four days 
assumed to have circular orbits. The longitudes of periastron and 
true anomaly of the first observation were taken to be uniformly 
distributed in the [0, 2 jt] interval and we adopted a uniform dis- 
tribution of the mass ratio q = ni2/mi between 0.1 and 1.0. The 
lower limit to the range of mass ratio was chosen because, in 
practice, it would be very difficult to observe the reflex motion 
of an O-star for such a low-mass companion. Finally, the orbital 
inclination was taken to be uniform in cos i between -1.0 and 
1.0. 

For four values of the primary mass (20, 30, 40, and 50 M Q , 
spanning roughly the range of spectral types for which some 
observational information on multiplicity is available) we sim- 
ulated a population of 100000 binary systems. For each system, 
we evaluated the maximum RV difference A RV = maxRV{<pi) - 
min RV(4>j) that would be measured on this system adopting the 
same temporal sampling as used for our five GIRAFFE-IFU ob- 
servation; *J The results are summarized in Tableland illustrated 
in Fig.fl4 If we assumed that a velocity difference of 20kms _1 
were required to claim the detection of a binary system^] we 
would find that fewer than 1.5% of the systems containing an 
04-6 V primary star (i.e. with a mass higher than 30 M Q ), and 
having an orbital period of shorter than 20 days, would escape 
detection. Only the lowest inclination systems would indeed re- 
main undetected. For longer orbital periods, the fraction of sys- 
tems that would be missed by our GIRAFFE-IFU campaign 
steeply increases: for systems with orbital periods between 50 

8 Our UVES campaign includes eight pointings. However, the three 
additional pointings only yield a marginal improvement in the detection 
efficiency compared to the results discussed here. 

9 None of the presumably single O-star objects discussed in this pa- 
per displays a velocity excursion of more than 20kms~' in the most 
stable photospheric lines. 



and 100 days, about half of the binaries would escape detection. 
However, owing to the assumed preponderance of short orbital- 
period systems, the total fraction of systems with a period up to 
100 days that we would miss, would still remain quite low (4 - 
5%). 
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Fig. 14. Distribution of ARV = max /? V(</>,) - minRV(^) for a 
population of massive binaries subdivided into various intervals 
of orbital period (labelled I, II, III, IV, and V, see Table[7]l. In this 
simulation, m\ was assumed to be equal to 40 M Q . 



The same approach was applied to the sampling of the three 
GIRAFFE-ARGUS observations (see Table|8]l. As could be ex- 
pected, this sampling is less efficient for the detection of bi- 
nary systems. Nonetheless, for primary masses in the range 20 
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- 50 M Q (corresponding roughly to 04-8 V spectral types), our 
simulations indicate that a binary system with a period shorter 
than 100 days has less than 10% a probability of escaping detec- 
tion because ARV < 20km s _1 . 

Therefore, we can safely conclude that the probability that 
either of the stars MSP 18, 171, 182, 183, 199, or 203 or stars A, 
B or C is a short -period binary system must be very low. If these 
stars are indeed binaries, they must have either a very low orbital 
inclination, a very long orbital period, or both. The situation is a 
bit less clear for WR 20b for which we do observe a maximum 
RV difference above our threshold. However, as stated above, 
this result must be considered with caution, since all the strong 
emission lines display line-profile variability that could easily 
be responsible for the observed shifts in RV, whilst the weaker 
emission lines are found to be more constant, in terms of both 
morphology and RV. 

The lack of a clear binary signature for WR20b and MSP 18 
is especially surprising given their high X-ray luminosities (see 
Naze et al. 120081) . These stars have logLx/^boi = -6.16 and 
-5.92 for WR20b and MSP 18 respectively, which are signifi- 
cantly larger than for typical single O-type stars and are as large 
as, or even larger than the value of the interacting wind system 
WR 20a (log Lx/Lboi = -6.15, Naze et al. 120081 ). Since colliding 
winds have the potential to generate copious amounts of X-rays 
(e.g. Stevens, Blondin & Pollock 1992), high X-ray luminosi- 
ties have traditionally been associated with interacting winds in 
a binary system. However, recent investigations have revealed 
that interacting wind binaries are not always X-ray overlumi- 
nous (e.g. Sana et al. 120061 Naze 120091 ), at least as far as O-star 
binaries are concerned, and that alternative mechanisms such as 
magnetically confined winds can also produce high X-ray lumi- 
nosities (e.g. Babel & Montmerle 1997 ). This implies that collid- 
ing winds, and by extension binarity, are neither a necessary nor 
a sufficient condition for high X-ray luminosities. Nonetheless, 
the high X-ray luminosities of MSP 18 and WR20b make these 
objects very interesting targets that deserve additional monitor- 
ing, in terms of both optical spectroscopy, to search for long- 
term RV variations, and spectropolarimetry, to search for evi- 
dence of magnetic fields or non-spherical winds. 

If these stars were indeed single, then their bolometric lumi- 
nosities would suggest masses in the range 50 - 70 M , provided 
that the main-sequence mass-luminosity relation holds for all 
these objects. It has been argued that the most massive stars 
might actually display an O If* supergiant or hydrogen-rich WN 
type spectrum rather than a normal O-type spectrum (Crowther 
et al. 120101 ). While it is true that the highest ever dynamical 
masses have been measured for WNLha stars so far (Rauw et 
al. EXH1 Niemela et al. l2008l Schnurr et al. 2008, 2009 ), it is 
puzzling that the WN6ha star WR20a and the 05V-III stars 
MSP 18 and MSP 203 have very similar bolometric luminosities 
and hence masses. In this context, it is worth recalling the 
result of De Becker et al. d2006l ). and Taylor et al. ( F20TT1 ). De 
Becker et al. (2006 ) derived a minimum dynamical mass of 
(150 + 50) M Q for th e 05.5 111(f) primary of HD 15558, whilst 
Taylor et al. (1201 II) inferred minimum dynamical masses of 
Mi sin 3 i = (78 ± 8) and M 2 sin 3 i = (66 + 7) M for the 
06.5 Iafc + 06Iaf binary system R 139. After all, these results 
suggest that some apparently 'normal' O-type stars could 
actually span a range of masses quite similar to that of WNLha 
stars. 

Although we have monitored only a limited sample of 
O-type stars in Westerlund2, our results suggest that unlike 
some other very young open clusters, such as NGC6231 (Sana 



et al. l2008l) or Trumpler 16 (Rauw et al. l2009l , Westerlund2 is 
most probably not very rich in short period early-type binaries. 
With only five binaries or binary candidates (including WR 20a 
and the three eclipsing OB systems, two of which host B-type 
primaries) out of the fifteen stars investigated, our results are 
more reminiscent of the situation encountered in IC 1805 (De 
Beck er et a l.[2006l Hillwig et al. [2006b and NGC 2244 (Mahy 
et al. 120091 ). If we restrict ourselves to O-type stars, the fraction 
of confirmed or candidate binaries drops to about 15% (2 out of 
13). Given the higher stellar density in the core of NGC 6231 
than in both NGC 2244 and IC1805, Mahy et al. d2009j ) 
tentatively suggested that the number of short-period early-type 
binary systems might be correlated with the cluster density. 
However, with its rather large cluster density and low fraction 
of short-period massive binaries, Westerlund2 apparently does 
not follow such a correlation. A more extensive multi-epoch 
monitoring of the other O-type stars in Westerlund2 would be 
highly valuable in helping us to establish the true binary fraction 
and investigate the question of a positional dependence of the 
occurrence of early-type binaries inside the cluster. 

The vast majority of the O-type stars in our sample indicate 
a RV of about 20 - 30km s as measured using the Hen 
A 5412 line. There are two exceptions to this rule, MSP 18 and 
MSP 171, which are instead found to have negative RVs for this 
line of (-1.1 ± 1.8)kms-' and (-9.3 ie^kms- 1 for MSP 18 
and MSP 171, respectively. Both stars are located outside the 
dense core of the cluster, at angular separations of 55 and 
23 "for MSP 18 and MSP 171, respectively. Therefore, these 
objects might have been ejected from the cluster core by means 
of dynamical interactions (see e.g. Leonard & Duncan [19901 ) or 
be long-period eccentric binaries observed near their quadrature 
phase. Alternatively, these stars could have strong stellar winds 
and their Hen A 54 12 line might then be blue-shifted owing to 
contamination by wind effects. This latter explanation is how- 
ever at odds with the absence of other strong wind signatures in 
their spectra (see the discussion of their spectra in the present 
work and in Paper I). We therefore conclude that these stars are 
runaway candidates or possible long-period eccentric binaries 
observed near a phase of maximum RV excursion. 

Finally, we have obtained the very first spectra of mid- to 
late-type O stars (MSP 223, A, B, C, D, and E) in Westerlund2. 
Together with the results of Paper I, the number of known O-type 
stars in the cluster is now 18. However, given the apparent mag- 
nitude of these mid to late-type O stars, only the tip of the ice- 
berg has yet been identified. Our current census of O-type stars 
in Westerlund2 is obviously still far from being complete and 
requires additional high spatial-resolution spectroscopic inves- 
tigations. Indeed, from the relationship between the mass of the 
most massive star and the mass of its parental cluster (Weidner et 
al. 120101 ), we estimate that if m max = 85 M B , Westerlund2 likely 
contains about 3000 M in the form of stars. Adopting the IMF 
of Weidner et al. (2010), we then estimate that the number of 
stars more massive than 20 M should be around 45. Therefore, 
more than half of the O-type stars of Westerlund 2 would still 
await discovery! 
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Appendix A: Tables of individual RV measurements 

The tables below provide the heliocentric RVs measured using 
various lines of stars that do not exhibit any clear binary signa- 
ture. 
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Table A.l. Radial velocities of WR20b as measured for the various emission lines in the spectrum. 



Date Radial velocity 



HJD-2454000 


H/J 


Nv/14945 


Nrvi5204 


Hen A 5412 


Nv A 5737 


Nrv/16383 


Hen A 6406 


Ha 


He ii ^6683 




(kms- 1 ) 


(kms- 1 ) 


(km s" 1 ) 


(kms" 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(km s- 1 ) 


(kms- 1 ) 


532.587 


132.2 


-22.0 


57.8 


108.0 


-52.9 


21.0 


-15.8 


71.1 


33.9 


532.631 


125.2 


-23.2 


33.6 


105.9 


-34.7 


-1.2 


-27.4 


71.5 


41.2 


534.581 


101.9 


-49.5 


36.6 


87.2 


-58.8 


28.9 


-30.2 


63.3 


3.1 


534.628 


91.2 


-20.9 


29.0 


78.7 


-48.3 


7.4 


-24.5 


68.9 


11.7 


536.624 


84.3 


-45.0 


49.0 


95.0 


-72.9 


46.8 


-15.3 


11.9 


44.8 


537.626 


125.7 


-39.4 


55.0 


105.1 


-38.5 


-14.2 


0.6 


43.5 


46.0 


538.719 


132.4 


-42.2 


74.3 


123.9 


-58.7 


13.5 


7.4 


100.5 


66.6 


539.622 


132.1 


-36.4 


64.8 


125.7 


-38.9 


46.4 


12.8 


125.8 


65.0 



Table A.2. Radial velocities of MSP 18 measured using various absorption lines and the weak C m A 5696 emission line. 



Date 






Radial velocity 






HJD-2454000 




Hen A 5412 


Cm A 5696 


Ha 


He ii A 6683 




(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


532.587 


-41.7 


-3.9 


-16.8 


-122.1 


4.2 


532.631 


-46.1 


1.7 


-0.9 


-123.8 


9.1 


534.581 


-45.5 


-2.1 


-17.0 


-122.3 


9.9 


534.628 


-33.5 


-2.0 


-17.3 


-118.4 


14.9 


536.624 


-42.8 


-2.0 


1.0 


-122.7 


6.3 


537.626 


-34.3 


0.6 


-8.3 


-120.3 


19.2 


538.719 


-50.1 


-0.5 


-6.8 


-122.3 


6.4 


539.622 


-44.6 


-0.3 


7.2 


-122.5 


0.8 



Table A.3. Radial velocities of MSP 171 measured using various absorption lines. 



Date Radial velocity 

HJD-2454000 Hen A 5412 O m A 5592 Civ ^5801 Civ ^5812 

(kms- 1 ) (kms- 1 ) (kms- 1 ) (kms -1 ) 

532.587 =Ol =405 ~T03 -13.5 

534.581 -14.4 -16.9 -7.8 -22.1 

536.624 -17.5 -2.6 -5.0 -17.4 

538.719 -8.6 -23.1 5.5 -21.7 

539.622 -5.8 -20.5 5.5 -18.3 



Table A.4. Same as Table A.3 but for MSP 182. 



Date 
HJD-2454000 



Hen A 5412 

(kms- 1 ) 



Radial velocity 
O mA 5592 C iv A 5801 



(kms- 1 ) 



(kms- 1 ) 



Civ ^5812 

(kms- 1 ) 



532.587 
534.581 
536.624 
538.719 
539.622 



34.5 
27.7 
33.3 
29.8 
27.0 



28.9 

30.4 
29.8 
35.6 



36.8 
36.5 
43.2 
38.0 
27.5 



4.1 
8.0 

1.3 
16.7 
-2.8 



Table A.5. Same as Table A.3 but for MSP 203 + MSP 444. 



Date 




Radial velocity 




HJD-2454000 


Hen A 5412 


Oiii,15592 


Civ A 5801 


Civi5812 




(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


532.587 


28.4 


47.7 


27.3 


-0.8 


534.581 


24.7 


38.8 


25.6 


5.5 


536.624 


23.3 


33.7 


21.1 


7.4 


538.719 


21.5 


39.4 


17.4 


0.0 


539.622 


23.9 


52.6 


19.8 


0.6 
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Table A.6. Radial velocities of MSP 183 measured using various absorption and emission (E) lines. 



Date 




Radial velocity 




HJD-2454000 


Hen ,15412 


Civ A 5801 (E) 


Civ /15812(E) 




(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


532.631 


23.0 


38.0 


53.9 


534.628 


25.3 


31.3 


55.5 


537.626 


23.5 


23.4 


49.1 



Table A.7. Radial velocities of MSP 199, Star B and C measured using the Hen A 5412 line. 



Date 


Radial velocity 




HJD-2454000 


MSP 199 Star B 


StarC 




(kms -1 ) (kms~') 


(kms -1 ) 


532.631 


29.1 18.0 


28.0 


534.628 


36.4 13.6 


29.7 


537.626 


26.2 21.3 


32.3 



Table A.8. Same as Table A.3 but for Object A. 



Date 






Radial velocity 






HJD-2454000 


He i A 5016 


Hen A 5412 


in A 5592 


Civ A 5801 


Civ A 5812 




(kms" 1 ) 


(kms" 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms" 1 ) 


532.631 


20.6 


26.3 


17.1 


23.0 


25.6 


534.628 


8.8 


28.1 


22.5 


22.0 


26.6 


537.626 


19.4 


28.5 


16.5 


18.7 


19.2 



